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Abstract 

Objectives: To evaluate the efficacy and safety of radiofrequency ablation (RFA) versus hepatic resection (HR) for early 
hepatocellular carcinoma (HCC) meeting the Milan criteria. 

Methods: A meta-analysis was conducted, and PubMed, Web of Science, the Cochrane Library, CBM, CNKI and VIP 
databases were systematically searched through November 2012 for randomized and nonrandomized controlled trials (RCTs 
and NRCTs). The Cochrane Collaboration's tool and modified MINORS score were applied to assess the quality of RCTs and 
NRCTs, respectively. The GRADE approach was employed to evaluate the strength of evidence. 

Results: Three RCTs and twenty-five NRCTs were included. Among 11,873 patients involved, 6,094 patients were treated 
with RFA, and 5,779 with HR. The pooled results of RCTs demonstrated no significant difference between groups for 1- and 
3-year overall survival (OS), recurrence-free survival (RFS) and disease-free survival (DFS) (p>0.05). The 5-year OS (Relative 
Risk, RR 0.72, 95% CI 0.60 to 0.88) and RFS (RR 0.56, 95% CI 0.40 to 0.78) were lower with RFA than with HR. The 3- and 5-year 
recurrences with RFA were higher than with HR (RR 1.48, 95% CI 1.14 to 1.94, and RR 1.52, 95% CI 1.18 to 1.97, respectively), 
but 1-year recurrence and in-hospital mortality showed no significant differences between groups (p>0.05). The 
complication rate (RR 0.18, 95% CI 0.06 to 0.53) was lower and hospital stays (Mean difference -8.77, 95% CI -10.36 to 
— 7.18) were shorter with RFA than with HR. The pooled results of NRCTs showed that the RFA group had lower 1-, 3- and 5- 
year OS, RFS and DFS, and higher recurrence than the HR group (p<0.05). But for patients with very early stage HCC, RFA 
was comparable to HR for OS and recurrence. 

Conclusion: The effectiveness of RFA is comparable to HR, with fewer complications but higher recurrence, especially for 
very early HCC patients. 
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Introduction 

Cancer is a major component of the global burden of disease 
(GBD). There were 2.49 billion disability-adjusted life years 
(DALYs), or 361 DALYs per 1000 population worldwide in 
20 1 0 [1] , and all neoplasms accounted for 7.6% (189 million 
DALYs) of global DALYs, an increase of 20 million DALYs 
(11.8%) compared with 2008[1,2]. A study based on the human 
development index (HDI) of Bray F et al [3] estimated an increase 
in the incidence of new cancer cases of all kinds to 22.2 million 
annually by 2030, and that increases would be proportionally 
greatest in low-HDI settings compared with high-HDI countries 
(76% vs. 25%). There were 19.1 million DALYs for liver cancer in 
2010, which accounted for 0.8% of the GBD or 10.1% of the 
DALYS for all neoplasms [1]. 



Hepatic resection (HR) and liver transplantation (LT) were 
recommended by the latest guidelines for early hepatocellular 
carcinoma meeting the Milan criteria, with the 5-year survival rate 
potentially reaching 50 to 75% [4,5]. However, only 20-35% of 
patients are suitable for liver resection because of the low diagnosis 
rate for early HCC and to poor liver function [6]. In addition, few 
patients can be treated with liver transplantation because of the 
strict inclusion criteria, high cost, and limited donor liver 
resources. 

Radiofrequency ablation (RFA) has higher efficacy and is 
associated with fewer complications and shorter hospital stays. 
RFA can also be administered repeatedly. Although RFA may 
gradually reach acceptability as an alternative treatment, the long- 
term efficacy and safety should still be evaluated systematically. 
Our previous study demonstrated that the overall quality of 
previously published systematic reviews and meta-analyses com- 
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paring RFA and HR for small hepatocellular carcinoma was poor, 
with an inadequate base of evidence [7] . Therefore, physicians may 
make an incorrect decision using these recommendations as best 
evidence without any additional quality evaluations to guide their 
clinical practice. 

The purpose of this study was to retrieve the best available 
evidence and produce a meta-analysis comparing the long-term 
results of RFA and HR for early hepatocellular carcinoma to 
reduce research bias and improve the quality of evidence. 

Materials and Methods 

Inclusion criteria 

The PICOS approach was used for eligibility criteria[8]: 

Population. Patients met the Milan criteria (single HCC 
smaller than 5 cm in diameter or up to 3 nodules that were each 
smaller than 3 cm in diameter) or the UCSF criteria (single tumor 
smaller than 6.5 cm in diameter or up to 3 nodules that were each 
smaller than 4.5 cm in diameter and 8 cm in total diameter) with 
liver function Child-Pugh class A or B (the number of patients with 
Child-Pugh C was no more than 10%). Patients were without 
major vascular invasion and lymphatic spread or extrahepatic 
metastasis, and had no previous treatment of HCC with any anti- 
cancer treatment [Transcatheter Arterial Chemoembolization 
(TACE), Percutaneous Ethanol Injection (PEI), and Microwave 
Ablation (MWA)]. 

Intervention. RFA 

Comparison. HR 

Outcome. Efficacy: O b overall survival rate (1-, 3-, 5-years), 
recurrence-free survival (RFS) (1-, 3-, 5-years), disease-free survival 
(DFS) (1-, 3-, 5-years). 

Safety: O z , mortality; recurrence rate (1-, 3-, 5-years); compli- 
cation rate 

Study design. Randomized controlled trials (RCTs), non- 
randomized controlled trials (NRCTs), retrospective clinical or 
cohort study 

DFS is defined as the time from randomization until recurrence 
of tumor or death form any cause [9]. In the case of HCC, the 
definition of DFS was identical as RFS where both recurrence and 
death form any causes are events [10]. However, the two terms 
were not clearly distinguished in the included RCTs or NRCTs, so 
we reported the two indicators respectively. 

Exclusion criteria 

Conference abstracts, reviews, letters, systematic reviews or case 
reports were excluded. Metastatic liver cancer (i.e., colorectal liver 
metastases) or recurrent hepatocellular carcinoma after resection 
was excluded. Those studies that mixed other effective interven- 
tions in either treatment group or control group (i.e., TACE, PEI, 
LITT) as well as any that had a length of follow-up of less than one 
year were also excluded. 

Data sources and search strategy 

We systematically searched 19 systematic reviews (SRs) 
comparing RFA with HR for small hepatocellular carcinoma in 
previously published studies and tracked the 39 primary studies 
included in these SRs. The six databases of PubMed, Web of 
Science, the Cochrane Library, CBM, CNKI, and VIP were 
systematically searched through November 2012. The following 
MeSH terms or free-words were used: ("hepatic resection [Title/ 
Abstract]" OR "surgical resection[Title/Abstract]" OR "liver 
resection [Title/ Abstract]" OR hepatectomy[MeSH Terms]) 
AND (radiofrequency [Title /Abstract] OR "radiofrequency abla- 
tion[Title/Abstract]" OR "catheter ablation[MeSH Terms]" OR 



"RFA[Title/Abstract]") AND ("hepatocellular carcinoma [Tide/ 
Abstract]" OR "liver neoplasm [MeSH Terms]" OR "liver 
tumor[Title/Abstract]" OR "liver cancer[Title/Abstract]"). We 
did not restrict the language of publication. 

Review selection and data extraction 

The PRISMA statement was followed when searching and 
screening the literature [8] . Two reviewers (WYQ^ LQOJ selected 
articles independendy by browsing titles and abstracts according to 
predefined inclusion and exclusion criteria. If necessary, judgment 
was made by viewing the full text. The two reviewers extracted 
data using standardized forms independently if the publication met 
the inclusion criteria. Discrepancies between the two reviewers 
were resolved by discussion or by a third person (LYP).The 
extracted contents included first author, publication year, type of 
study, inclusion criteria, number of participants, age, gender, 
number of nodules, tumor size, the Child-Pugh score, length of 
follow-up, and clinical outcomes. 

Quality assessment 

Two reviewers (WYQ LQOJ evaluated the methodological 
quality of included studies independendy. The tool for evaluating 
the risk of bias in the Cochrane Handbook for Systematic Reviews 
of Interventions [1 1] and the modified MINORS scores question- 
naire^] were applied to assess the quality of RCTs and NRCTs, 
respectively. GRADE profiler 3.6 was employed to evaluate the 
strength of the evidence. 

Statistical analysis 

Meta-analysis was conducted using RevMan 5.1 software. For 
dichotomous variables, the relative risk (RR) and the odds ratio 
(OR) with 95% confidence intervals (CIs) were used for RCTs and 
NRCTs, respectively, For continuous variables, mean difference 
(MD) with 95% CI was applied. The fixed-effect model (Mantel- 
Haenszel) was used if the result of the heterogeneity test was 
p>0.05 and if I 2 <50%. Otherwise, the random-effect model was 
applied. A p value <0.05 was considered significant. 

Subgroup analysis 

Subgroup analysis was applied based on predetermined 
subgroup factors (i.e., tumor size, number of nodules, and Child- 
Pugh class). 

Sensitivity analysis 

We excluded the studies that might cause heterogeneity and re- 
estimated the combined effect of values, and then compared the 
results with the primary outcomes to verify the robustness of the 
outcomes. 

Publication bias 

The funnel plot and Egger's test were applied by Revman 5. 1 
and Stata 10.0 software, respectively. If the result of the test was 
p<0.05, it suggested that potential publication bias may exist. 
Otherwise, publication bias was considered absent. 

Results 

Search results 

After initial screening, 76 studies were identified. Of these, 48 
studies were removed after viewing full-texts for various reasons: a) 
mixed with other effective interventions (i.e., TACE, PEI) in either 
intervention group or control group (14 articles), b) without control 
group (11 articles), c) recurrent HCC (8 articles), d) did not meet 
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3S7 of records identified through database 

starching: 

.Pubmed(n»$6) 

.Web of science (rpl6) 

. Cochrane librarv(n«5) 

.CB^n-77) 

.CNKXn-156) 

.YTP(n-17) 



39 of additicnal records identified through 
pervious Nfeta-analysis, Systematic re\ie\v or 
Health Technolcev Assessment searched before 
201! 



352 of records afier duplicates removed 



T 



352 of records screened 



276 of records excluded after checking titles and 
abstracts 



■m 



76 of &ll-text articles assessed 
for elieibilttv 



2S of studies included for quality 
assessment and meta-analysis 
.3 RCT 
-25NRCT 



4S of full-text articles excluded, with reasons; 

. 14 Accepted other effective interventions 

.11 No control group 

.SRecuirent HCC 

. 7 Not meets the purpose of studies 

.4 Lack of detail baseline information 

. 3 Liver metastases 

. 1 Abstract only 



Figure 1. PRISMA flowchart of searching and selecting guidelines. 

doi:1 0.1 371 /journal. pone.0084484.g001 



the purpose of this study for other reasons (7 articles), e) lack of 
detailed baseline information (4 articles), f) metastatic liver cancer 
(3 articles), and g) abstract only (1 article). Finally, 28 studies 
including 3 RCTs [10,13,14] and 25 NRCTs [6,15-38] published 
between 2004 and 2012 were included (see Figure 1, Figure 2). 



Baseline characteristics of included studies 

A total of 1 1,873 patients with primary HCC were included in 
the 28 studies. Among those, 8,567 cases (72%) were males, and 
6,094 cases were treated with RFA and 5,779 cases with hepatic 
resection. There were 11,251 cases (94.8%) with Child- Pugh class 



■= 5 





Figure 2. Bibliometric map of included studies. 

doi:10.1371/journal.pone.0084484.g002 
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A, 1,632 cases (13.7%) with Child-Pugh class B, and only 17 cases 
(0.1%) with Child-Pugh class C. Patients had an average age of 50 
years, with varying degrees of HBsAg-positivity and cirrhosis. The 
mean length of follow-up ranged from 10 to 60 months (Table 1). 

Blinding was not performed in the three RCTs, and the overall 
quality of evidence level was grade B [10,13,14]. The quality of the 
NRCTs was moderate, with an estimated mean MINORS score 
(18 of total) of 15.8 (95% CI, 15.4-16.3). Only 15 (60%) studies 
were scored >16 [6,15-38] (Table 1). 

Clinical outcome of HCC patients with tumor size smaller 
than 5 cm 

Overall Survival. The pooled meta-analysis from the three 
RCTs [10,13,14] showed no significant difference of 1- and 3-year 
survival rates between groups (RR 0.98, 95% CI: 0.89-1.09, 
NNH = 33.3, p = 0.71; and RR 0.98, 95% CI: 0.74-1.29, 
NNH = 22.2, p = 0.87, respectively) (level of evidence: moderate). 
Only Huang et al[10] demonstrated that the 5-year survival rate in 
the RFA group was lower than in the HR group (RR 0.75, 95% 
CI: 0.60-0.88, NNH = 4.8, p = 0.001) (level of evidence: high) 
(Figure 3, Table 2). 

The pooled meta-analysis from the NRCTs showed that the 1-, 
3- and 5-year survival rates in the RFA group were significandy 
lower than in the HR group (OR 0.78, 95% CI: 0.63-0.97, 
NNH= 166.7; OR 0.67, 95% CI: 0.52-0.85, NNH= 12.5; and 
OR 0.58, 95% CI: 0.36-0.94, NNH= 10.6, respectively) (level of 
evidence: very low) (Table 3). 

Recurrence-free survival. The meta-analysis from the 
RCTs showed no significant difference of 1- and 3-year 
recurrence-free survival rates between groups (level of evidence: 
moderate to high). The 5-year survival rate in the RFA group, 
however, was lower than in the HR group (RR 0.56, 95% CI: 0.40 
to 0.78, NNH = 4.4) (level of evidence: high) (Figure 4, Table 2). 

The pooled results from the NRCTs demonstrated that the 1-, 
3- and 5-year recurrence-free survival rates in the RFA group were 
lower than in the HR group (OR 0.78, 95% CI: 0.64-0.95, 
NNH = 25.6; OR 0.67, 95% CI: 0.56-0.79, NNH = 8.5; and OR 
0.63 95% CI: 0.40-1.00, NNH = 7.5, respectively) (level of 
evidence: very low). However, there was no significant difference 
between groups of 5-year recurrence-free survival rates for patients 
with Child-Pugh class A (OR 0.64, 95% CI: 0.35-1.17; 
NNH = 8.3) (level of evidence: very low) (Table 3). 

Disease-free survival. Only one RCT [14] reported the 
disease-free survival rate, showing no significant difference 
between groups (level of evidence: moderate). 

The pooled meta-analysis from the NRCTs showed that the 1 -, 
3- and 5-year disease-free survival rates in the RFA group were 
significandy lower than in the HR group (OR 0.46, 95% CI: 0.38- 
0.55, NNH = 7.2; OR 0.49, 95% CI: 0.34-0.69, NNH = 5.1; and 
OR 0.52, 95% CI: 0.32-0.84, NNH = 5.9, respectively) (level of 
evidence: very low to low) (Table 3). 

Recurrence. The pooled results from the RCTs showed that 
the 3- and 5-year recurrence rates in the RFA group were 
significandy higher than in the HR group (RR 1.48, 95% CI: 
1.14-1.94, NNH = 7.1; and RR 1.52, 95% CI: 1.18-1.97, 
NNH = 4.6, respectively) (level of evidence: moderate). However, 
there were no significant differences of 1-year recurrence rates or 
of local or distant recurrence rates between groups (level of 
evidence: moderate) (Figure 5, Table 2). 

The results of meta-analysis for the NRCTs showed that the 1-, 
3- and 5-year recurrence rates in the RFA group were higher than 
in the HR group (OR 1.5, 95% CI: 1.03-2.19, NNH= 12.4; OR 
1.87, 95% CI: 1.23-2.84, NNH = 6.4; and OR 2.34, 95% CI: 
1.76-3.1 1, NNH = 5.9, respectively) (level of evidence: very low to 
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Study or Subgroup 



RFA 
Events Total 



HR 

Events Total 



Risk Ratio 
Weight M H. Random, 95'. CI 



Risk Ratio 
M-H, Random, 95*. CI 



1 -year survival rate 
Chen MS 2006 
Feng K 201 2 
Huang JW 2010 
Subtotal (95S CI) 
Total events 
Heterogeneity: Tau= = 0.01 ; Chi* = 11.48. df 
Test for overall effect: Z = 0.37 (P = 0.71) 



68 
81 
1 00 

249 



n 

84 
115 
270 



84 
78 
113 



90 
84 
115 
289 



33.3% 
33.6% 
33.1% 
100.0S 



1 .03 [0.95, 1 .1 0J 

1.04 [0.97. 1.12] 
0.88 [0.82, 0.95) 
0.98 [0.89,1.09] 



275 



2 (P= 0.003); P =83% 



3-year survival rate 



Chen MS 2006 


58 


71 


66 


BO 


33.3% 


1.11 [0.94.1.32] 


Feng K 201 2 


63 


84 


56 


84 


31.9% 


1.13(0.93.1.37] 


Huang JW 2010 


30 


115 


106 


115 


34.8% 


0.75 [0.66, 0.86] 


Subtotal (95S CI) 




270 




289 


100.0 : : 


0.98 [0.74.1.29] 


Total events 


201 




238 









Heterogeneity: Tau 2 = 0.05; ChP= 17.93, df= 2 (P= 0.0001); P= 89% 
Test for overall effect: Z = 0.1 7 (P = 0.87) 



5-year survival rate 

Huang JW 2010 63 115 87 115 100.0% 0.72 [0.60, 0.88] 

Subtotal (95% CI) 115 115 100.0% 0.72(0.60.0.88] 

Total events 63 87 

Heterogeneity: Not applicable 
Test for overall effect: Z= 3.23 (P = 0.001) 

1 — I — I 1 — I 

0.5 0.7 1 1.5 2 
Favours HR Favours RFA 



t 



Figure 3. Overall survival (OS) at 1-, 3- and 5-year in RCTs. 

doi:1 0.1 371 /journal.pone.0084484.g003 



low). However, there were no significant differences between 
groups for patients with Child-Pugh class A (level of evidence: very 
low) (Table 3). 

In-hospital mortality, complication rate and length of 
hospital stay. The pooled results of the RCTs showed no 
significant differences of in-hospital mortality between groups, but 
the complication rate in the RFA group was lower than in the HR 
group (RR 0.18, 95% CI: 0.06-0.53, NNT = 3.5) (level of 
evidence: moderate). Length of hospital stay in the RFA group 
was 8.77 days fewer than in the HR group (95% CI: 10.36 to 7.18 
lower) (level of evidence: moderate) (Figure 6, Table 2). 

The pooled results of the NRCTs for in-hospital mortality and 
complication rate were similar to those of the RCTs, but the 
hospital length of stay in the RFA group was 6.74 days fewer than 
in the HR group (95% CI: 1 1.33 to 2. 14 lower) (level of evidence: 
very low) (Table 3). 

Clinical outcome of HCC patients with tumor size 
between 3 and 5 cm 

The pooled results of the RCTs showed no significant difference 
of overall survival rates between the groups for single HCC 
patients with tumor size ranging from 3 to 5 cm in diameter (level 
of evidence: low to moderate) (Table 4). 

The results from the NRCTs showed that the 5-year survival 
rate in the RFA group was lower than in the HR group (OR 0.43, 
95% CI: 0.25-0.73, NNH = 5.3), but there were no significant 
differences of 1-and 3-year survival rates between groups (level of 
evidence: very low). The 1-, 3- and 5-year disease-free survival 
rates in the RFA group were lower than in the HR group (OR 
0.47, 95% CI: 0.26-0.83, NNH = 6.9; OR 0.35, 95% CI: 0.18- 
0.67, NNH = 6.1; and OR 0.18, 95% CI: 0.05-0.61, NNH = 9.9, 
respectively) (level of evidence: very low to low) (Table 5). 



Clinical outcome of HCC patients with tumor size smaller 
than 3 cm 

The pooled results of RCTs showed no significant differences of 
the overall survival rates between groups for patients with a single 
HCC and tumor size S3 cm in diameter (level of evidence: low to 
moderate) (Table 4). 

The pooled meta-analysis from NRCTs showed no significant 
differences of 1- and 3-year survival rates between groups, but the 
5-year survival rate in the RFA group was lower than in the HR 
group (OR 0.62, 95% CI: 0.43-0.90, NNH = 4.5; p = 0.01) (level 
of evidence: very low). The 1-, 3- and 5-year disease-free survival 
rates in the RFA group were lower than in the HR group, but the 
5-year disease-free survival rates for patients with a single HCC 
showed no significant difference between groups (OR 0.69, 95% 
CI: 0.35-1.36). The 1-, 3- and 5-year recurrence rates also showed 
no significant differences between groups, but the 5-year 
recurrence rate in the RFA group was lower than in the HR 
group for patients with Child-Pugh class A (OR 0.42, 95% CI: 
0.19-0.93, NNT = 5.1) (level of evidence: moderate). The compli- 
cation rate in the RFA group was lower than in the HR group 
(OR 0.44, 95% CI: 0.29-0.65, NNT=11.2) (level of evidence: 
low), but there was no significant difference between groups for 
patients with Child-Pugh class A (Table 6). 

Clinical outcome of solitary HCC patients with tumor size 
smaller than 2 cm 

The pooled results of three NRCTs [15, 16,22] showed that there 
was no significant difference in 1-, 3- and 5-year overall survival or 
RFS, 3-, 5-year DFS, and 1-,3-year recurrence between 
groups(p>0.05). However, the 1-year DFS (OR 0.22, 
95%CI:0.07-0.65, NNH = 4.4; p = 0.06), 5-year recurrence(OR 
0.42, 95%CI:0. 19-0.93, NNT = 5.0;P = 0.03) and complication 
rate (OR 0.23, 95%CI:0. 1 1-0.49, NNT = 3.1; p = 0.0001) in RFA 
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Study or Subgroup 



RFA 
Events Total 



HR 

Events Total 



Risk Ratio 
y.v-:nht ;.--■<. Random. 95 : : CI 



Risk Ratio 
M-H. Random, 95% CI 



1-year recurrence-free survival 

Feng K 201 2 76 84 72 84 51.1% 

Huang JW 2010 94 115 98 115 48.9% 

Subtotal (95* CI) 199 199 100.0% 

Total events 1 70 1 70 
Heterogeneity: Tau 3 = 0.00; Chi'= 1 .41, df= 1 (P= 0.24); P= 29% 
Test for overall effect: Z= 0.1 5 (P = 0.88) 



1.06 [0.94, 1.18] 
0.96(0.85.1.08] 
1.01 [0.92,1.11] 



3-year recurrence-free survival 

Feng K 201 2 51 84 42 84 49.1% 1.21 [0.92, 1.60] 

Huang JW 2010 53 115 70 115 50.9% 0.76 [0.59, 0.97] 

Subtotal (95% CI) 199 199 100.0% 0.95 [0.60, 1.52] 

Total events 1 04 112 

Heterogeneity: Tau' = 0.09; Chi'= 6.31, df= 1 (P= 0.01); l'= 84% 

Test for overall effect: Z = 0.20 (P= 0.85) 



5-year recurrence-free survival 

Huang JW 2010 33 115 59 115 100.0% 0.56 [0.40, 0.78] 

Subtotal (95% CI) 115 115 100.0% 0.56 [0.40, 0.78] 

Total events 33 59 

Heterogeneity: Not applicable 

Test for overall effect: Z= 3.36 (P = 0.0008) 



H 1 1 1 h 

0.2 0.5 1 2 5 
Favours HR Favours RFA 



Figure 4. Recurrence-free survival rate (RFS) at 1-, 3- and 5-year in RCTs. 

doi:10.1371/journal.pone.0084484.g004 



group were lower than that of HR group(P<0.05)(level of 
evidence: very low)(Table 7). 

Sensitivity analysis 

Among the three RCTs, the patients included by Huang et 
al[10] were the oldest (average age of 56 years). With exclusion of 
this study, the 1- and 3-year survival rates still showed no 



significant difference between groups, but the results of the 
heterogeneity test (I 2 ) changed from 83-89% to 0%. The study of 
Feng et al[13] included more patients with cirrhosis and a high 
proportion (61.5% to 75%) of patients with more than two 
nodules. With exclusion of this study, the complication rate in the 
RFA group was still lower than in the HR group (RR 0.12, 95% 



RFA HR Risk Ratio Risk Ratio 

Study or Subgroup Events Total Events Total Weight M-H, Fixed, 95% CI M-H. Fixed. 95% CI 



1-year recurrence 




12 84 


8 


84 


36.4% 


1 .50 [0.65, 3.48] 






Feng K 201 2 






1 — ■ 


Huang JW 2010 




19 115 


14 


115 


63.6% 


1.36(0.72, 2.57] 






Subtotal (95% CI) 




199 




199 


100.0% 


1.41 [0.85, 2.35] 






Total events 




31 


22 












Heterogeneity: Chi" = 


0.03, df=1 (P = 0.85); I s = 


0% 










Test for overall effect: Z = 


1.32 (P = 0.19) 














3-year recurrence 


















Feng K 201 2 




42 84 


32 


84 


55.2% 


1.31 [0.93,1.86] 






Huang JW 2010 




44 115 


26 


115 


44.8% 


1.69 [1.12,2.55] 




■ 


Subtotal (95% CI) 




199 




199 


100.0% 


1.48(1.14,1.94] 






Total events 




86 


58 












Heterogeneity: Chi" = 


0.88, df=1 (P = 0.35); l 5 = 


0% 










Test for overall effect: Z = 


2.90 (P = 0.004) 














5-year recurrence 


















Huang JW 2010 




73 115 


48 


115 


100.0% 


1.52(1.18.1.97] 




t 


Subtotal (95% CI) 




115 




115 


100.0% 


1.52 [1.18,1.97] 






Total events 




73 


48 












Heterogeneity: Not applicable 














Test for overall effect: Z = 


3.20 (P = 0.001) 










t 1 1 


I 1 1 



0.2 0.5 1 2 5 
Favours RFA Favours HR 



Figure 5. Recurrence rate at 1-, 3- and 5-year in RCTs. 

doi:10.1371/journal.pone.0084484.g005 
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Study or Subgroup 



Experimental 
Events Total 



Control 
Events Total 



Risk Ratio 
Weight M-H, Random, 95*. CI 



Risk Ratio 
M-H, Random. 95*. CI 



Chen MS 2006 
Feng K 2012 
Huang JW 2010 



71 
84 
115 



50 
18 
32 



90 
84 
115 



29.9% 
36.3% 
33.9% 



Total (95S CI) 270 289 100.0*. 

Total events 16 1 00 

Heterogeneity: Tau"= 0.64; Chi'= 7.76, 01= 2 (P = 0.02); P= 74% 
Test for overall effect: Z = 3.1 5 (P = 0.002) 



0.08 [0.02, 0.23] 
0.44 [0.20. 0.97] 
0.16 [0.06. 0.39] 

0.18 [0.06, 0.53) 



0.001 0.1 1 10 1000 
Favours RFA Favours HR 



Figure 6. Complication rate of RCTs. 

doi:10.1371/journal.pone.0084484.g006 



CI: 0.06-0.24), but the results of the heterogeneity test (I 2 ) 
changed form 74% to 0%. 

Publication bias 

As is shown in Table 8, we only analyzed the publication bias 
for indicators included in 10 or more studies[l 1]. After viewing the 
funnel plot and Egger's test, it was found that the indicators of 5- 
year survival for patients with tumor size smaller than 5 cm or 
3 cm in diameter showed no publication bias (P>0.05). The 
remaining indicators all showed that some degree of publication 
bias existed. 

Discussion 

HR is considered the preferred treatment for patients meeting 
Milan criteria with a single nodule or multiple lesions with good 
liver function but unsuitable for liver transplantation [4,5]. One 
prospective RCT demonstrated that the 3-year survival rate with 
liver resection for small hepatocellular carcinoma could reach 
74.8% [13]. However, 80% of cases were unsuitable for liver 
resection for various reasons such as low rate of early diagnosis, 
poor expected function of residual liver after surgery, and 
anticipated serious post-operative complications [39], Local 
ablation with RFA or PEI is recommended by the latest updated 
EASL-EORTC guidelines as the standard care for patients with 
Barcelona clinic liver cancer (BCLC) 0 to A level but were 
unsuitable for surgery in 20 12 [4]. 

RCTs of high quality and with a large sample size are 
considered the best sources of evidence, but RCTs are accompa- 
nied by high costs, are difficult to conduct, and often have poor 
external validity. RCTs performed in the field of surgery, 
especially with double-blinding methods, are extremely difficult. 
Systematic review and/or meta-analysis of RCTs are often the 
most feasible methods to address this situation. Previous meta- 
analyses worldwide have compared the long-term efficacy of RFA 
and HR for the treatment of small HCC, but prospective RCTs 
with large sample sizes were rarely included. The systematic 
review authors usually combined the results of RCTs and NRCTs 
together or mistook retrospective studies as RCTs, leading them to 
draw conclusions of low reliability [39-43], as various risks for bias 
or confounding factors might have existed among the study 
designs. However, the impact these biases on outcomes could not 
be ascertained because of the limited information provided in 
these primary studiesfl 1]. Therefore, such results without a more 
strict risk assessment would inevitably lead to an erroneous 
estimation of effects and mislead clinical decision-making. 

In this study, three RCTs without blinding were included, all 
with only moderate quality of evidence [10,13,14]. The pooled 
results of our meta-analysis showed no significant differences 
between the RFA and HR groups in overall survival and 
recurrence-free survival rates at 1 and 3 years and in recurrence 



rates at 1 year following the treatment of small hepatocellular 
carcinomas meeting the Milan criteria. The RFA group had 
higher recurrence rates at 3 and 5 years and lower complication 
rates when compared with the HR group. It is well-known that 
tumor size, number of lesions, location, liver function, the presence 
of portal vein invasion, the presence of vascular invasion and the 
width of the tumor-free margin during surgical excision were the 
independent prognostic factors affecting the survival of pa- 
tients [44,45]. The recurrence rate after RFA is related to 
incomplete tumor ablation. In case of HCC in which local 
curative ablative therapy was obtained by securing a safety 
margin [46]. It is considered that the local recurrence rate 
markedly differed among patients with or without a sufficient 
safety margin [47]. However, RFA is a minimally invasive 
procedure, and it is often hard to achieve a specific safety margin 
in three dimensions all around a large tumor [46], leading to a 
higher recurrence but lessened compromise of liver function, 
which might be one of the reasons for the relatively low 
postoperative complication rate. Commonly, patients with small 
HCC rarely dies within 1 year and the recurrence impacts the 
overall survival gradually. However, the 1-year survival rate of 
RFA in the RCTs Huang et al[10] conducted was 87%, whereas 
that of resection was 98 % . This result is markedly different from 
the other two RCTs [13,14], which might be the factors leading to 
the clinical heterogeneity. There are some reasons may contribute 
to this problem. Firstiy, of 108 RFA-treated patients in the study of 
Huang et al[10], seventeen patients had lesions in dangerous 
locations, leading difficult to secure a safety margin. Secondly, the 
average age of included patients in both groups were much old, 
and they were prone to the therapy of resection. Apart form that, 
the rate of loss to follow-up was greater in the resection group (18/ 
115, 15.6%) than in RFA group (7/115, 6.1%), which would 
probably influence the comparison between two groups. 

NRCTs of high quality with large sample sizes could provide 
research evidence for a wider population, particularly with a 
greater advantage for assessing the safety of an intervention. 
However, NRCTs might also be more easily influenced by 
different kinds of bias or by unknown confounding factors [1 1]. In 
this study, twenty-five NRCTs with an average quality of 
"moderate" and a large pooled sample size of 11,314 subjects 
were included. The pooled results of the meta-analysis could have 
some significance for guiding clinical practice as they showed that 
the overall survival in the RFA group was significantly lower than 
in the HR group for patients with small tumor size less than 5 cm 
in diameter. However, identical results were suggested otherwise in 
respect of recurrence, complication rate, in-hospital mortality and 
hospital length of stay between the RCTs and NRCTs. The reason 
for this phenomenon was that rare RCTs with low heterogeneity 
and large sample size were included which showed no difference 
between the groups, but did not rule out the impact of potential 
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Table 8. Publication 


bias of studies which were 


more than ten. 










Indicators 


Subgroup 


Years 


No. of studies 


Egger's test P value 


Publication bias 


Overall survival 


HCCS5 cm 


1-y 


23 


-C0.001 


Yes 






3-y 


23 


<0.001 


Yes 






5-y 


15 


0.182 


No 




HCCS5 cm and Child A 


1-y 


11 


0.001 


Yes 






3-y 


12 


0.006 


Yes 




HCCS3 cm 


1-y 


14 


0.003 


Yes 






3-y 


15 


0.005 


Yes 






5-y 


12 


0.21 


No 


Disease-free survival 


HCCS5 cm 


1-y 


12 


0.039 


Yes 






3-y 


11 


0.011 


Yes 


Complication rate 






15 


<0.001 


Yes 



doi:1 0.1 371 /journal.pone.0084484.t008 



confounding factors or bias (i.e., publication bias) in the NRCTs, 
leading to a possible overestimation of the survival effect in the 
NRCTs. 

Few studies except for three retrospective ones have reported 
the outcome of solitary HCC with tumor size less than 2 cm. 
Nevertheless, the conclusion of these studies was varies. WangJH 
et al[15] and Hung HH et al [2 2] have concluded that RFA was as 
effective as HR in patients with BCLC very early stage HCC, 
while Peng ZW et al[16] showed that percutaneous RFA was 
better than those of HR, especially for central HCC. In this study, 
we found that the overall survival and 1-, 3-year recurrence in 
RFA group were as equal effective as HR group, and the 
complication rate and 5-year recurrence in RFA group were lower 
than HR group. However, the overall quality of evidence 
evaluated by GRADE (the Grading of Recommendations 
Assessment, Development and Evaluation) was very low to low 
due to limited sample size or many known or unknown risk of bias 
existed in observational studies. Therefore, when applying the 
outcome of the evidence to clinical practice, physicians should take 
caution. More large-scale, well-conducted RCTs or retrospective 
studies focused on the topics were still needed in the further. 

GRADE (2011 version) has identified five categories or factors 
downgrading of the quality of evidence: study limitations (risk of 
bias), imprecision, inconsistency, indirectness, and publication bias 
[48,49]. For observational studies, the GRADE group has also 
identified three categories or factors raising the quality of evidence: 
large magnitude of effect, dose-response gradient, and plausible 
confounding, all of which can increase confidence in the estimated 
effects [50]. Therefore, our study strictly followed the quality 
assessment criteria of the GRADE guidelines, evaluated the 
quality of evidence for the important outcomes of each patient, 
and performed a subgroup analysis based on risk factors (i.e., the 
Child-Pugh class, tumor size, and number of nodules) that were 
likely to affect the clinical outcome, so as to reduce the impact of 
risk factors on clinical outcomes to some extent. It is well-known 
that cirrhosis is also an independent prognostic factor for patients 
with HCC. For example, in the study of Nishikawa H et al [20], 
they concluded that the presence of liver cirrhosis was the sole 
significant factor for recurrence-free survival. However, based on 
the current primary studies, it is difficult to apply a subgroup 
analysis for this factor individually. We suggest that more high 
quality RCTs or retrospective studies should be focus on this topic 
in the further. 



Limitations 

Only three RCTs were included, and the limited sample size 
and high heterogeneity might affect the robustness of the clinical 
outcomes to some extent. Given the characteristics of the included 
NRCTs, we deleted some items (i.e., prospective collection of data, 
unbiased and blinded assessment of the study endpoint, and 
prospective calculation of the study size) from the MINORS score 
questionnaire which were unsuitable for this study, but we were 
convinced that it would not have an impact on the quality of 
evidence of the final clinical outcome. 

Conclusions 

The pooled meta-analysis of the RCTs demonstrated no 
significant difference between groups of the overall survival, 
recurrence-free survival, disease-free survival and in-hospital 
mortality rates for early HCC with tumor size smaller than 
5 cm in diameter, but the RFA group had higher recurrence rates, 
lower complication rates, and shorter hospital lengths of stay. 

The pooled meta-analysis of the NRCTs showed no significant 
difference in recurrence rates between groups for patients with 
Child-Pugh class A or tumor size smaller than 3 cm or 3 cm to 
5 cm in diameter. However, when combining the results of all 
patients with tumors smaller than 5 cm in diameter, it is showed 
that the RFA group had lower overall survival and higher 
recurrence rates. For patients with Child-Pugh A/B and tumor 
size less than 2 cm in diameter, the pooled results concluded that 
RFA was as effective as HR in the overall survival and 1 -, 3-year 
recurrence, and RFA yielded lower complication rate and 5-year 
recurrence than HR. However, it is still need RCTs of high quality 
to further enhance the level of evidence. 

Based on full consideration of the current available best 
evidence, a comprehensive conclusion can be drawn: RFA is 
comparable to HR with lower complication rates but with higher 
recurrence rates. All relevant risk factors that may affect the final 
outcome of patients should be considered, so as to balance 
minimizing recurrent HCC after RFA with improving the quality 
of life of patients. The authors suggest that more high quality 
RCTs or retrospective studies should focus on RFA versus HR for 
very early HCC patients and to provide the better clinical decision 
for physicians. 
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